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Economy is demanding new models, able to understand and predict the evolution of 
markets. To this respect, Econophysics offers models of markets as complex systems, that 
try to comprehend macro-, system- wide states of the economy from the interaction of 
many agents at micro-level. One of these models is the gas- like model for trading markets. 
This tries to predict money distributions in closed economies and quite simply, obtains 
the ones observed in real economies. However, it reveals technical hitches to explain 
the power law distribution, observed in individuals with high incomes. In this work, non 
linear dynamics is introduced in the gas-like model in way that an effort to overcome these 
flaws. A particular chaotic dynamics is used to break the pairing symmetry of agents 
0,0- The results demonstrate that a "chaotic gas-like model" can reproduce 
the Exponential and Power law distributions observed in real economies. Moreover, it 
controls the transition between them. This may give some insight of the micro-level 
causes that originate unfair distributions of money in a global society. Ultimately, the 
chaotic model makes obvious the inherent instability of asymmetric scenarios, where 
sinks of wealth appear and doom the market to extreme inequality. 

Keywords: Econophysics; Money Dynamics; Chaotic Systems; Complex Systems; Com- 
putational models 

PACS Nos.: 11.25.Hf, 123.1K 

1. INTRODUCTION 

After the economic crisis of 2008, much criticism has been thrown upon modern 
economic theories. They have failed to predict the financial crisis and foresee its 
depthless [[T] . On the whole, this failure has been charged to the inherent complexity 
of the markets [ [2] . 

In this respect, Econophysics has been ultimately offering new tools and per- 
spectives to deal with economic complexity [ [3] - [5] . Simple stochastic models have 
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been developed, where many agents interact at micro level giving rise to global em- 
pirical regularities observed in real markets. These models are giving some guidance 
to uncover the underlying rules of real economy. 

One of the most relevant examples of these models is the conjecture of a kinetic 
theory with (ideal) gas-like behavior for trading markets [|4]-|9]. This model offers 
a simple scheme to predict money distributions in a closed economic community of 
individuals. In it, each agent is identified as a gas molecule that interacts randomly 
with others, trading in elastic or money-conservative collisions. Randomness is also 
an essential ingredient of this model, for agents interact in pairs chosen at random 
and exchange a random quantity of money. In the end, the model shows that the 
asymptotic distribution of money in the community will follow the exponential 
(Boltzmann-Gibbs) law for a wide variety of trading rules [[5]. 

This result may elucidate real economic data, for nowadays it is well established 
that income and wealth distributions show one phase of exponential profile that 
covers about 90 — 95% of individuals (low and medium incomes) [[10]- H2]- Despite of 
this fact, the model shows some technical hitches to explain the other phase observed 
in real economies. This is the power law (Pareto) distribution profile integrated by 
the individuals with high incomes [[12]- [15]. The model needs to introduce additional 
elements to obtain this distribution, such as saving propensity among traders or 
diffusion theory [ [5] , [9] . 

The work presented here, intends to contribute in this respect. To do that, it 
proposes to incorporate chaotic dynamics in the traditional gas-like models. This is 
done upon two facts that seem particularly relevant in this purpose. The first one is 
that, determinism and unpredictability, the two essential components of chaotic sys- 
tems, take part in the evolution of Economy and Financial Markets. On one hand, 
there is some evidence of markets being not purely random, for most economic 
transactions are driven by some specific interest (of profit) between the interacting 
parts. On the other hand, real economy shows periodically its unpredictable com- 
ponent with recurrent crisis. The prediction of the future situation of an economic 
system resembles somehow the weather prediction. Therefore, it can be sustained 
that market dynamics appear to be chaotic; in the short-time they evolve under de- 
terministic forces though, in the long term, these kind of systems show an inherent 
instability. 

The second fact is that the transition from the Boltzmann-Gibbs to the Pareto 
distribution may require the introduction of some kind of inhomogeneity that breaks 
the random indistinguishability among the individuals in the market. Though this 
is not a necessary condition, it may produce it. By now there have been differ- 
ent approaches that obtain this kind of transition, in [ 16 a random market with 
inhomogeneous saving propensity obtains the Pareto distribution. The use of de- 
terministic dynamics for modeling of financial markets has also been undertaken 
[ [17] and deterministic dynamics with homogeneous models are also able to obtain 
the Pareto distribution [ 18 . In this work, a deterministic model with asymmetric 
chaotic transactions is considered. The dynamics that governs the market is going 
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to be selected so that is able of breaking the pairing symmetry of interacting agents 
(J, i). A chaotic map with parametric-controlled symmetry may be an ideal 
candidate to obtain that, as it can quite simply produce parametric-controlled rou- 
tines of evolution for the market. Additionally, a chaotic system will also represent 
a simple mechanism to introduce two other important ingredients: some degree of 
correlation between agents and establishing a complex pattern of transactions. 

These concepts or hypothesis have inspired the introduction of a model for trad- 
ing markets with chaotic patterns of evolution. Amazingly, it will be seen that a 
chaotic market is able to reproduce the two characteristic phases observed in real 
wealth distributions. The aim of this work is to observe what may happen at a micro 
level in the market, which is responsible of producing these two global phases. 

This paper is organized as follows: section 2 introduces the basic theory of the 
gas-like model and describes the simulation scenario used in the computer sim- 
ulations. Section 3 shows the results obtained in these simulations. Final section 
gathers the main conclusions and remarks obtained from this work. 

2. SIMULATION SCENARIO OF THE CHAOTIC MARKET 

The model proposed here is a multi-agent gas-like scenario [ H] . The study of these 
scenarios was first proposed by the authors in [ [19] . There, it is shown that the use 
of chaotic numbers produces the exponential as well as other wealth distributions 
depending on how they are injected to the system. This paper considers the sce- 
nario where the selection of agents is chaotic, while the money exchanged at each 
interaction is a random quantity. 

As in the gas-like model scenario a community of N agents is given an initial 
equal quantity of money, too- The total amount of money, M = N*mo, is conserved. 
The system evolves for a total time of T = 2 * N 2 transactions to reach the asymp- 
totic equilibrium. For each transaction, at a given instant i, a pair of agents 
with money (m^rrij 1 ) is selected chaotically and a random amount of money A to 
is traded between them. The amount Am is obtained through equations (1) where 
v is a float number in the interval [0, 1] produced by a standard random number 
generator: 

Am = v * (toj* + m/)/2, 
TO; t+1 = m,* — Am, 
mj t+1 = mj l + Am. 

This particular rule of trade is selected for simplicity and extensive use. Com- 
parisons can be established with popularly referenced literature [H]-[9]. Here, the 
transaction of money is quite asymmetric as agent j wins the money that i losses. 
Also, if agent i has not enough money (to/ < Ato), no transfer takes place. 
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Fig. 1. (a) Representation of 2000 points of the chaotic attractor of the Logistic bimap. This 
bimap is symmetric respective to the diagonal y=x when A a = Aj, . (b) Representation of the 
spectrum of coordinate xt .The values of the parameters for both graphics are A a = A = 1.032 in 
the blue curves and X a = 1.032, A;, = 1.084290 in the black curves. The curves show the symmetric 
case, A a = A5 = 1.032 and the most asymmetric one, while chaotic behaviour is still observed, 
A a = 1.032, A 6 = 1.084290 

The chaotic selection of agents for each interaction, is obtained by a par- 
ticular 2D chaotic system, the Logistic bimap, described by Lopez-Ruiz and Perez- 
Garci'a in [ 20 . This system is given by the following equations and depicted in 

Figru 

T : [0, 1] x [0, 1] — > [0, 1] x [0, 1] 

xt=\ a {3yt-i + l)xt-i(l-x t -i), (2) 
y t = \b(3x t -i + l)y t _i(l - t/t-i). 

The dynamical properties of this system make it an ideal candidate for our 
proposes. The model presented here requires a chaotic system able to break the 
pairing symmetry of interacting agents As it is explained below, the 

symmetry of this chaotic map can be controlled parametrically and so, this will 
allow to produce series of symmetric or asymmetric agents interactions at will. 

The chaotic system of equation [2] considers two Logistic maps that evolve in a 
coupled way. This system is used to select the interacting agents at a given time. 
The pair is easily obtained from the coordinates of a chaotic point at instant 
t, X t = [xt,yt], by a simple float to integer conversion: 

i = (int)(x t * N), 

(3) 

j = (int)(y t * N). 

This procedure will make that the selection of a winner ( j ) , or a looser (i) , follows 
a chaotic pattern. Moreover, as there are two Logistic maps are coupled with each 
other, this bimap is able to introduce a strong correlation in the selection of agents 
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of each group. This correlation is modulated or governed by the parameters A. The 
Logistic bimap presents a chaotic attractor in the interval X a ,b G [1.032,1.0843] 
(see a real-time animation in [[21]). Consequently, the chaotic selection of agents 
is guaranteed by taking some appropriate A a ,b in this range. FigQ] (a) shows two 
trajectories for two different groups of values of A a and Ah, showing maximum 
asymmetry while prevailing chaotic behaviour. 

A property of the Logistic bimap is its symmetry respective to the diagonal 
y = x , when A a = A& (see Fig[T] (a)). The spectrum of coordinate Xt shows a peak 
for w = 0.5 (see Fig[T](b)) presenting an oscillation of period two that makes a jump 
over the diagonal axis alternatively between consecutive instants of time. Both sub- 
spaces x > y and y > x are visited with the same frequency and the shape of the 
attractor is symmetric. When A& becomes greater than A a the part of the attractor 
in sub-space x < y becomes wider and the frequency of visits of each sub-space 
becomes different. These characteristics are used in the model as an input variable, 
so the degree of symmetry in the selection of agents can be modified at will. 

In the end, the dynamics of the economic interactions are going to be complex 
and quite different from the random scenario, where any agent may interact with 
any other with the same probability. Here, market transactions are restricted in 
the sense that one agent will only interact with specific groups of other agents. To 
see this, just think that an i agent given by an x — xq coordinate, with xo any 
constant value in the interval [0, 1]. Draw the line x — xq in FigfT] (a) and then it 
will be seen that this agent will interact with its neighbors on the diagonal (j i) and 
maybe with two other distant groups of agents. Additionally these secondary groups 
interact with other groups and so on, giving rise to a complex flow of interactions in 
the whole market. This may seem more realistic than the random case, as normally 
one individual doesn't interact directly with any other, but with specific groups of 
traders, being the complex connections of the community, the ones that develop a 
global indirect trading. 

3. SIMULATION RESULTS OBTAINED IN THE CHAOTIC 
GAS-LIKE MODEL 

Taking into account the chaotic market model discussed in the previous section, 
different computer simulations are carried out. In these simulations a community of 
N = 5000 agents with initial money of mo = 1000$ is taken. The simulations take 
a total time of T — 2 * N 2 = 50 millions of transactions. 

Different cases are produced when different values of the chaotic parameters A a 
and A;,. In this way the symmetry of the selection of agents can be varied at will. 
This will allow us to see the effects of introducing an asymmetry in the process 
of selection of interacting agents. Simulations are produced when A a = 1.032 has a 
fixed value and Ab varies from 1.032 to 1.033162 with increments of S = 0.581 xlO -4 . 

The results of these simulations reveal very interesting features. First, there is 
a group of individuals that keep their initial money and don't interact at all. This 
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Fig. 2. Final CDF's obtained in the simulations for different axis scales, (a) Final CDF's for 
parameters A a = 1.032 and X b = 1.032, 1.324648, 1.326391, 1.329296 and 1.331039 in natural units, 
(b) Representation of the same CDF's in a logarithmic plot, (c) Representation of the same CDF's 
in a double logarithmic plot. 

result would be consistent with real markets where not every agent are active. This 
number is due to the shape of the chaotic attractor. In Fig[T](a), it can be seen that 
the attractor hardly reaches the extreme values of the interval [0,1]. In the case 
of A a = Xb = 1.032 the number of non interacting agents is 1133 in a community 
of N = 5000. When Afc increases the attractor expands and this number becomes 
smaller. 

After removing passive agents and their money, the final distributions of money 
are shown in Fig. [2] with different axes scales and for different values of A;,. Here 
A a = 1.032 for all cases. As it can be seen, the symmetric case produces an expo- 
nential distribution, as in the random scenario of the gas-like model [5] . Amazingly, 
increasing the asymmetry of the chaotic selection, the money distribution degener- 
ates progressively from the exponential shape to a power law profile. 

Figj2] (a) shows in five curves, the cumulative distribution functions 
(CDF's) obtained for the following four cases: A a = 1.032 and A;, = 
1.032, 1.324648, 1.326391, 1.329296 and 1.331039. Here the probability of having a 
quantity of money bigger or equal to the variable MONEY, is depicted in natural 
axis plot. The symmetric case A a = Ah = 1.032 is highlighted in blue color. Fig|2] 
(b) shows the same data in a natural logarithmic plot, to illustrate the exponential 
profile obtained in the symmetric case (highlighted in blue). It also shows how a 
progressive asymmetry in the selection of agents degrades this characteristic profile. 
Fig|2] (c) shows again the same data but in a double-logarithmic plot. It gives an 
extensive view of the degradation observed before, making evident the straight line 
power-law profile obtained in the most asymmetric case (highlighted in blue) . 

As a consequence, one may say that a chaotic selection of agents is reproducing 
the two characteristic distributions observed in real wealth distributions. On one 
hand, one obtains the exponential distribution, as in the random gas-like models. 
This characteristic appears when chaotic market is operating under a symmetric rule 
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of selection of interacting agents. The significance may be that the correlations in 
agent interactions looks like random, and all agents win or lose with fair probability. 

On the contrary, when asymmetry is introduced in the chaotic process of selec- 
tion of agents, the distribution of money becomes progressively more unequal. The 
probability of finding an agent in the state of poorness increases and only a minority 
of agents reach very high fortunes. What is happening here is, that the asymme- 
try of the chaotic map is selecting a set of agents preferably as winners for each 
transaction (j agents). While others, with less chaotic luck, become preferably the 
losers (i agents) . As in real life, there are markets where some individuals possess a 
preferred status and this makes them win in the majority of their transactions. 

From these primary results, it seems interesting to study the dynamics of the 
system at micro detail. This might help to uncover the possible causes that originate 
unfair distributions of money in society. Then, simulations are repeated for a more 
tractable number of agents N = 500. The initial money is mo — 1000$ as previously, 
and the simulations take a total time of T = 2 * N 2 = 0.5 millions of transactions. 

First the CDF is obtained for three simulation cases with different A& values 
(see Fig[3](a)). The degradation of the exponential distribution is again obtained 
as asymmetry increases. However when Fig[2] (a) and Figj3] (a) are compared an 
avalanche effect becomes now evident for bigger markets. Here one may appreciate 
that the shape of the money distributions in the symmetric case is the same for 
a simulation scenario of N — 5000 or N = 500 agents. In contrast, the asymmet- 
ric cases produce more unequal or dramatic distributions when the community of 
agents is greater. A slight change in the asymmetry of the interactions, compare 
for example the case A& = 1.033162, produces a significantly higher degradation 
in the final distribution of money for the N = 5000 market. At this point, it is 
interesting to remark that, in a chaotic scenario the size of the market matters. A 
small asymmetry in the rules of selection of agents will become far more aggressive 
in its final distribution of wealth depending on the size of the market. 

As a consequence, the globalization of markets that evolve under chaotic rules, 
may suffer avalanche effects that drive more unequal and dramatic distributions of 
wealth. This may resemble some situations observed in real economy, where a small 
unbalance in global markets, may produce large differences in the share of capital 
among individuals. 

Secondly, the network of chaotic economic interactions is analyzed. Agents trade 
in a complex network of interactions illustrated in Fig|3](b). As we can see, agents 
in the range 320 to 360 have quite a higher number of interactions between them. 
Amazingly, these agents are not specially treated in the symmetric case, for in this 
scenario all agents have the same opportunities. On the other hand, the figure shows 
how in the most asymmetric case some agents are preferably selected as winners. 
The maximum asymmetry case depicted here treats some agents as winners in all 
transactions. See highlighted in yellow, agent 30, who is never selected as an i agent, 
a looser. In this situation the group of the agents in the range 320 to 360, which 
interact mainly among themselves, will always get poorer as the majority of the 
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Fig. 3. Results obtained for a simulation scenario with N = 500 agents, (a) Representation of 
the same CDF's of Fig[2]but for N = 500 agents and three different cases (X b = 1.032, 1.033162 
and 1.08429). (b) Number of transactions between pairs of agents for the most asymmetric case 
(A a = 1.032 and A;, = 1.08429). Axis x and y represent pairs of agents selected as winner 

(j) or loser (i) respectively. Axis z represents the number of transactions of each pair. Agent 30 
position is highlighted with a yellow line to show that he never looses (i agent), (c) Time evolution 
of money for agent 348 in the symmetric case (A a = A;, = 1.032), and agents 348 and 30 in the 
most asymmetric case (A;, = 1.08429). 

society is getting poor. The more they interact the sooner they will lose all their 
money. Put it in another way, in the asymmetric case, not all agents have the same 
opportunities. 

In third place, FigJ3] (c) shows the evolution in time of an agent's money in 
symmetric and asymmetric cases. Here agent number 348, who in the range 320 to 
360, is depicted as an example. In the picture it can be observed that when the 
market is symmetric, agents' wealth oscillates in a random-like way. Becoming rich 
or poor in the end it is just a question of having a very profitable transaction, being 
in the right place at the right time. In the chaotic market model this is simply the 
result of some specific selection of the initial conditions and chaotic parameters. Let 
us note that in the symmetric case, the money distribution becomes exponential. 
On a global perspective, there are no sources or sinks of wealth and even then, when 
the individuals have the same opportunities, the final distribution is unequal. 

However when the market is asymmetric agent 348 becomes inevitably poorer 
as time passes. In this case there are agents that never lose (as for example, agent 
30) and they become inevitably sinks of wealth. This can be observed in Fig[3] (c) 
where agent 30 is precisely one of these sinks, accumulating more and more money 
every time he trades. These chaotic favored agents make the rest of the community 
poorer. On a global perspective, becoming rich in the asymmetric case depends on 
who you are and who you interact with (the shape of the chaotic attractor and 
its initial conditions). This resembles a society where some individuals belong to 
specific circles of economic power. 

Finally, to illustrate these qualitative observations, we study the number of times 
that agents win or lose in the simulation cases. The results obtained are depicted in 
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Fig. 4. Representation of the role of N = 500 agents in the community after all interactions. 
Agents are arranged in descending order according to their final wealth (agent is the richest 
of all). The upper graphic shows the total number of wins over looses of an agent. The bottom 
graphic shows the number of times an agent has been selected as i agent or looser, (a) Symmetric 
case with X a = A5 = 1.032. Here, the number of wins and looses is uniformly distributed among 
the community. A range of 112 individuals are passive and always have the initial 1000$. (b) 
Asymmetric case with A a = 1.032 and A5 = 1.084290. Here, 46 agents never interact and 18 agents 
never loose. 



Fig]?] Here it is shown the number of times (interactions) that an agent has been a 
looser (bottom graph) and the difference of winning over losing times (top graph). 
The x axis shows the ranking of agents arranged in descending order according to 
their final wealth. So, agent number is the richest of the community and agent 
number 500 is the poorest. 

FigH](a) depicts the symmetric case, where A a = \ = 1.032. Here, the number 
of wins and looses is uniformly distributed among the community. There is also a 
range of agents that don't interact (112 agents), this can be seen clearly in the figure 
now. In this case, the chaotic selection of agents shows no particular preference for 
any other and the final distribution becomes the exponential. This is similar to 
traditional gas- like simulations with random agents [ [4] - [9] . 

FigH] (b) shows the same magnitudes for the asymmetric case with A a = 1.032 
and Af, = 1.0842908. Here the asymmetry is maximum prevailing chaotic behavior. 
In this case, there is a group of agents in the range of maximum richness that never 
loose. The chaotic selection is giving them maximum luck and this makes them 
richer at every transaction. A range of agents, lower than in the symmetric case, 
are passive and never interact (46 agents). In this case, the final money distribution 
becomes quite unequal: 269 agents (half of the society approximately) end half 
poorer with a final wealth inferior to 500$ and of them, and 56 agents (approx. the 
10%) finish with less that 50$. It is also interesting to see in Fig|4] (a) that in the 
poor class, there are agents that have a positive difference of wins over looses, but 
amazingly they become poor anyway. Consequently, one can deduce that they are 
also bad luck individuals. They may be selected as j agents in most part of their 
transactions, but unfortunately their corresponding trading partners (i agents) are 
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poor too, and they can effectively earn low or no money in these interactions. 

To summarize, these results show that the symmetric case of selection of agents 
resembles a society where agents have equal opportunities in the market. In this 
case, agents are equally selected winners or losers (symmetry, x < — > y and i < — > j) 
in a chaotic and correlated fashion. This resembles a market where all agents have 
the same opportunities of success. In this case, there is no particular preference 
for any other and the final distribution of money becomes exponential, as in the 
random gas-like model. 

Quite the opposite, when a small asymmetry is introduced in the rule of selecting 
agents, the society becomes more unequal. Some agents become preferred winners 
for most of their transactions and the money. The opportunities of winning in the 
end just depend on interacting with the proper group of rich agents. This resembles 
a market where there are established circles of power and corruption. Additionally, 
the chaotic nature of the market may also reproduce an avalanche effect as a result of 
globalization of markets, where a small asymmetry may produce extreme inequality. 

4. CONCLUSIONS 

This work introduces a novel approach in the field of economic (ideal) gas-like 
models [H]-|S]. It proposes a chaotic selection of agents in a trading community and 
it accomplishes to obtain the two most important money distributions observed in 
real economies, exponential and power-law distributions, besides a mechanism of 
transition between them. 

The introduction of chaos is based on two considerations. One is that real econ- 
omy shows a kind of a chaotic character. The other one is that dynamical systems 
offer a very simple and flexible tool able to break the inherent symmetry of random 
processes and overcome the technical difficulties of gas models to shift from the 
Boltzmann-Gibbs to the Pareto distributions. 

A specific 2D dynamical system under chaotic regime is considered to produce 
the chaotic selection of agents. This system represents two logistic maps coupled 
in a multiplicative way and produces a strong correlation between the interacting 
agents. This correlation guaranties a complex behavior that can be tuned through 
two system parameter ( A a and A&) and so, it is able to produce symmetric or 
asymmetric conditions in the market. 

The results obtained in this model show how a chaotic selection of agents under 
symmetric conditions produces a final exponential distribution of money. In con- 
trast, when asymmetry is introduced in the selection of agents, the distribution of 
money becomes progressively more unequal until it produces power law profiles. 

The analysis of micro level interactions in both cases shows what is happen- 
ing in the market and what is responsible of producing these two global phases 
(exponential or power law). 

The main conclusions produced by this model may resemble some situations in 
real economy. In a symmetric scenario the selection of agents resembles a society 
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where agents have equal opportunities in the market. In this case, there is no par- 
ticular preference for any other and the final distribution becomes the exponential. 
This is similar to traditional gas-like simulations with random agents. On a global 
perspective, there are no sources or sinks of wealth and even when the individuals 
have the same opportunities, the final distribution is unequal. From an agent point 
of view, becoming rich or poor just depends of being at the right time in the right 
place. 

In the asymmetric scenarios, a small asymmetry is introduced in the rule of se- 
lecting agents and then, the society becomes unequal. Some agents become preferred 
winners for most of their transactions. The opportunities of winning in the end just 
depend on being a preferred agent and interacting with the proper group of rich 
agents. The final money distribution resembles the Pareto distribution. Moreover, 
the chaotic nature of the market also reproduces the impact of the size of the market 
on the final distribution of wealth. When unequal or asymmetric conditions of trad- 
ing are ruling, an avalanche effect may be produced and then inequality becomes 
extreme. 

The authors hope that this work that may bring new ideas and perspectives 
to Economy and Econophysics. The proposal of considering chaotic dynamics in 
multi-agents modelling may also be of interest to other fields, where scientists try 
to describe and understand complex systems. 
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